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Abstract. - We find the distribution of transmission eigenvalues in a series of identical junc- 
tions between chaotic cavities using the circuit theory of mesoscopic transport. This distribution 
rapidly approaches the diffusive wire limit as the number of junctions increases, independent 
of the specific scattering properties of a single junction. The cumulant generating function and 
the first three cumulants of the charge transfer through the system are obtained both in the 
normal and in the superconducting state. 



Introduction. - Higher-order correlators of current fluctuations in mesoscopic conduc- 
tors have been studied extensively over the last decade both theoretically [1-7] and experi- 
mentally [8-18]. The reason is that they contain, in general, additional information to the 
usual differential conductance such as higher moments of the transmission eigenvalue dis- 
tribution, the value of effective charge involved in transport processes, the size of internal 
energy scales of the system or the correlations intrinsic to the many-body state of entangled 
systems [19-21]. While the conductance is proportional to the average transmission probabil- 
ity of the structure at low temperatures, the current noise power Pi depends on the second 
moment of transmission eigenvalue distribution which is characterized by the Fano factor 
F = Pj/2el — [J2 n T n (l — T„)]/ J2 n Tn- Here e is the electron charge, I is the average current 
through the sample, and T„ are the transmission eigenvalues. Recent experiments on noise 
confirmed the theoretical predictions [2,3] on the universal distributions of transmission eigen- 
values in a metallic diffusive wire [8,9] and in an open chaotic cavity [10] with Fano factors 
F = 1/3 and F — 1/4, respectively. The crossover from a single cavity to the diffusive wire 
limit as the number of internal junctions increases was studied experimentally by Oberholzcr 
et al. [10] and Song et al. [11] recently. 

Particle-hole correlations introduced by a superconducting terminal also modify the noise. 
The low temperature noise of the subgap transport is doubled for tunnel junctions [12] and 
in diffusive normal wires in contact with a superconductor [13, 14]. The noise in an open 
cavity is found to be more than two times larger in the superconducting state [15] than in the 
corresponding normal state junction, in agreement with theoretical predictions [2]. 
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The third correlator contains the first three moments of transmission eigenvalue distribu- 
tion and is related to the asymmetry of the distribution [22] . In contrast to the current noise 
which is thermally dominated at temperatures larger than the bias voltage according to the 
fluctuation-dissipation theorem, the third correlator is in this regime proportional to the cur- 
rent, without the need to correct for the thermal noise. However, higher-order correlators are 
increasingly more difficult to measure because of the statistical fluctuations [23] and the influ- 
ence of environment [17,20]. Recent measurements of the third-order correlations of voltage 
fluctuations across the nonsuperconducting tunnel junctions by Bomze et al. [18] confirmed 
the Poisson statistics of electron transfer at negligible coupling of the system to environment. 

The statistical theory of transport, full counting statistics [24,25], provides the most de- 
tailed description of charge transfer in mesoscopic conductors. The semiclassical cascade 
approach to higher-order cumulants based on the Boltzmann-Langevin equations has been 
developed by Nagaev et al. [26]. The stochastic path integral theory of full counting statis- 
tics was introduced by Pilgram et al. [27]. The quantum-mechanical theory of full counting 
statistics based on the extended Keldysh-Green's function technique [28-30] in the discretized 
form of the circuit theory [31,32] was put forward for multiterminal circuits by Nazarov and 
Bagrets [33]. In this article we use the circuit theory [34,35] to study the elastic quasiparticle 
transport in arrays of chaotic cavities focusing on the crossover from a single cavity to the 
universal limit of a diffusive wire [36] as the number of inner contacts increases. We find the 
analytical expressions for the distribution of transmission eigenvalues, the cumulant generating 
function and the first three cumulants both in the normal and in the superconducting state, 
generalizing the previous results on noise [10] in such a system to all higher-order correlations. 
The similar finite-size effects on the noise and the third correlator have been studied numer- 
ically by Roche and Dougot [37] within an exclusion model. Ballistic to diffusive crossover 
in metallic conductors with obstacles as a function of increasing disorder has been studied 
by Macedo [38] within the scaling theory of transport combined with the circuit theory. The 
effects of Coulomb interaction on the current and noise in chaotic cavities and diffusive wires 
have been studied by Golubev et al. [39]. The noise in series of junctions has been measured 
by Oberholzer et al. [10] and Song et al. [11] recently. 

Transport in an array of chaotic cavities. - The system we consider consists of chaotic 
cavities in series between TV identical junctions characterized by N c h transverse channels with 
transmission eigenvalues {T n }. We can neglect the energy dependence of the transmission 
eigenvalues of the system if the electron dwell time is small with respect to time scales set by 
the inverse temperature and applied voltage. Also we neglect the charging effects assuming 
that the conductances of the junctions are much larger than the conductance quantum 2e 2 /h. 
The quasiparticle distribution function is isotropic between junctions due to chaotic scattering 
in the cavities. We apply the circuit theory of mesoscopic transport and represent the specific 
parts of the system by the corresponding discrete circuit elements, as shown in fig. ^ The 
Greens functions of the leads are denoted as Go(0) = Gl(0) and Gn(x) = Gr(x), while the 
Greens functions of the internal nodes are Gi(x), i = 1, . . . , N — 1. The counting field x can 
be incorporated through the boundary condition [35] at the right lead according to 

G N ( X ) = e -' l{x/2)fK Gjv(0) e i{x/2)fK . (1) 

Here Go,n(0) are the bare Greens functions of the Fermi leads in the Keldysh(~) <g) Nambu(~) 
space, ?i and 5", are the Pauli matrices and fx = Ti&$- The connection between adjacent 
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Fig. 1 - An array of chaotic cavities in series between the Fermi leads (top) and the discrete circuit- 
theory representation of the system (bottom). The leads and the cavities are associated with the 
corresponding matrix Greens functions Gi. The junctions are characterized by a set of transmission 
eigenvalues {T„}. The matrix current is conserved throughout the circuit. 



nodes is described in general by a matrix current [31] 



2e^ 
h 



N ch 

E 



2T„[G.;+i, Gi 



A + T n ({G i+1 ,Gi}-2y 



(2) 



which flows from the node i to the node i + 1. The set of circuit-theory equations for the 
Greens functions of the internal nodes consists of matrix current conservations Ii,i+i = const 
and normalization conditions Gf = 1. As shown in ref. [40], we can seek for the solution in the 
form Gi — FiL({Go, Gn}) Gq + Fm({Go, Gn}) Gn, where the functions FiL,n depend only 
on the anticommutator {Go, Gn} and therefore commute with all Gj. As a consequence, the 
anticommutators {Gi, Gj} depend only on {Go, Gat} and commute with all G&. We emphasize 
that the above consideration is independent of the concrete matrix structure of the Greens 
functions and relies only on the quasiclassical normalization conditions G\ = 1. Because the 
junctions are identical, the matrix current conservation reduces to [40] 



Gi = 



Gi 



G^ 



{Gi-i, Gi+i} + 2 



(3) 



Taking the anticommutator of eq. © with Gi and Gj+i , respectively, we find that {Gj_i, Gi} — 
{Gi, Gi+i} = Q' , for all i. Our aim is to find Q' in terms of Q = {Go, G/v}. Now we take the 
anticommutator of eq. © with Go and obtain the following difference equation 



7i+i - Q'li + Ji-i = 0, 



(4) 



where ji = {Go, GT. After solving eq. (J2J with the boundary conditions 70 = 2 and = G, 
and using that 71 — G' , we find 



G' = {(G + VG 2 - 4)/2]V^ + [(g - VS 2 - 4)/2] 1 / ff . 



(5) 



The cumulant generating function S(x) of charge transfer through the structure can be 
obtained as a sum of the actions of the connected pairs of nodes [33]. For identical junctions 
in series S(x) reduces to the contribution of a single junction multiplied by N: 

S(x) = ~ J d£ tr[5( X )], with S( X ) = N ^ ln[l + T n <$' - 2)/4]. (6) 
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Fig. 2 - Distribution of transmission eigenvalues for TV open contacts in series (left panel) and for 
contacts of lower transparency T n — 0.4 and T„ < 1 (right panel), normalized to the transmission 
distribution of a diffusive wire pd(T). Ballistic to diffusive crossover in a metallic disordered wire is 
shown for comparison (left panel, dash-dotted curves) as a function of increasing disorder L/l = TV — 1. 



Here to is the total measurement time which is much larger than the characteristic time scale 
on which the current fluctuations are correlated. For a large number of junctions, TV 3> 1, 
the cumulant generating function given by eq. © approaches the universal limit of a diffusive 
wire [25,32] S(x) = (—tog/8h) J d£ tr[arccosh 2 (Cy/2)], which does not depend on the specific 
scattering properties {T n } of a single junction, the shape of the conductor or the impurity 
distribution [34,36]. Here g = (J2 n T n )/N is the total conductance of a wire in units of 2e 2 /h. 

The distribution of transmission eigenvalues of a composite junction £>jv.{t„}C7") is directly 
related to the cumulant generating function by [32] 

1 ^fdS(C) 



PN,{T n} (T) = ^ Im(— - J, (7) 



d( C=-l/T-i0 

where £ = (Q — 2)/4. From eqs. ©, 10 and Q we find 

, sin(7r/iV) 

PN,{T n }(J) = PD(T) — jjy 

4T»[(fr 1 / jY + r 1 / jv )(2 - T n ) + 2T w cos(7r/AQ] 

[(fci/iv + 6-i/JV)( 2 - T n ) + 2T n cos(7r//V)] 2 - 4(1 - T^b 1 ^ - f 



(8) 



Here pd{T) = (,g/2)(l/TVl — T) is the transmission distribution of a diffusive wire, b = 
(1 + Vl — T) 2 /T, and (• • • ) = (J2 n ■ ■ ■ )/Q>2 n T n ) denotes the averaging over the transmission 
eigenvalues of a single junction. For N open contacts in series (T n = 1) eq. JSJ reduces to 

PN( ' PD{ ' tt/N (l + % /T3T)2/A r + (l- VT^T)2/iv + 2Ti/^cos(7r//V)' (> 

The crossover from a single cavity to the diffusive regime as the number of junctions N 
increases is shown in fig.Elfor open contacts (left panel) and for contacts of lower transparency 
T n = 0.4 and T n <^ 1 (right panel). The transmission distribution of a metallic disordered 
wire of the length L and the mean free path I is shown for comparison (dash-dotted curves in 
ng.E cf. ref. [3]). 

The transmission eigenvalue distributions given by eqs. J3J and @ can be probed ex- 
perimentally by measuring higher-order correlators of current fluctuations across the junc- 
tion at low temperatures. The first three moments of charge transport statistics are re- 
lated to the average current, the current noise power and the third correlator according to 
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I = i(e/t )d x S\ x = , Pi = (2e 2 /t )d x S\ x=0 , and d = -i(e 3 /t )^S'| x= o, respectively. In 
the linear regime, which we consider here, the current is proportional to the bias voltage 
with conductance given by g — E„T„)/JV = f Q dTp N ^ Tn j(T) T in units of 2e 2 /h. At 
temperatures much lower than the bias voltage, the current noise power and the third corre- 
lator are linear in the current, with the slopes given by the Fano factor F = dPi /d(2el) = 
(1/g) /„* dTp N , {Tn} (T)T(l-T) and the "skewness" C = dd/die 2 !) = (1/g) J* dTp N>{Tn} {T) 
x T(l — T)(l— 27"), respectively. For the normal-state junction the two parameters are given 



3 V N 2 

^ r _ 1 , 5(2 -3(T 2 )) 4-30(T„ 2 (l-r„)) \ 

C_ 15V. A 2 + A 4 " J' (H) 

The Fano factor given by eq. (|10|) coincides with the result previously obtained within Boltz- 
mann-Langevin approach which takes into account both cavity noise and partition noise at 
the contacts and was confirmed experimentally for up to three open contacts in series [10]. 
The sign of C is related to the asymmetry of transmission distribution [22], being negative 
(positive) when the weight of the distribution is shifted towards open (closed) transmission 
channels. Equation (| 1 II) shows that closed channels prevail in the composite junction for 
N > 2 even for completely open inner contacts, in agreement with eq. 

Now we focus on the junction sandwiched between a normal-metal and a superconductor in 
the coherent regime in which we can neglect the particle-hole dephasing (E t h ^> |eV| , fcsT, A 
where E t h is the inverse dwell time). At temperatures and bias voltages smaller than the 
superconducting gap A, the transport properties can be obtained by integrating the An- 
dreev reflection probability Ra = T 2 /(2 — T) 2 over the transmission distribution [2,40] and 
correcting for the effective charge e* = 2e. The conductance, the Fano factor and the skew- 
ness are given by gs = dT pN,{T n }(X) Ra (in units of Ae 2 /h), Fs — dPi/d(2e*I) = 
(l/g s ) Jo dTp N>{Tn} (T)R A (l-R A ), and C s = dd/d(e* 2 I) = (l/g s ) ft dTp N>{Tn} (T)R A (l- 
R A )(1 — 2Ra), respectively. For N open contacts in series we find 



N ch 1 _ 1 

9S = 7777 ¥7 /.atM F S=-Z 



~ -2 



and 



CS 15 



2A cos 2 (tt/4A)' 3L 4 A 2 Vcos 2 (tt/4A) 

1 



(12) 



1 



2 ) + 8Nl{ 2 + 15 co S HnHN) )\- (13) 



4 A 2 Vcos 2 (7r/4A) I 8A 4 V cos 4 (tt/4A) J. 

The skewness C$ given by eq. \1'6\ is positive for A > 1 which indicates that the asymme- 
try of the Andreev reflection distribution Pna{Ra) = Pn{T)(IT / 'cLRa is in favor of closed 
Andreev channels [40]. For the more general case of N contacts characterized by transmis- 
sion eigenvalues {T n } we obtain g s = (J2 n at n )/N and Fs = 1 - (£ n "nAiVC a ™)' wnere 
a„ = VRjiv^Rj + cos(vr/2A)]/[l + V^cos(tt/2A)] 2 , /3„ = 2/3 - (1/6A 2 ){1 - 6a„ + 
Sv^f/IV^f + cos(vr/2A)]}, and R% = T 2 /(2 - T„) 2 . Fano factors and skewnesses for the 
normal-state and superconducting junctions are shown in fig. as a function of the number 
of contacts in series and for different contact transparencies. It is interesting to note that 
in the coherent superconducting regime, which we consider here, the higher-order correlators 
satisfy the approximate scaling relations F$(N) w F(2N) and Cs(N) ~ C(2N) which are 
exact for incoherent Andreev transport [41]. For large A this results in the full reentrance of 
transport properties of a diffusive wire in contact with a superconductor [30] as a function of 
the particle-hole coherence. 



6 



EUROPHYSICS LETTERS 




2 lj 1 1 1 1 li 1 1 1 1 J 0.2 0.0 — ■ ' 1 1 1 1 1 1 J 0.0 

23456234567 23456234567 

N N N N 



Fig. 3 - The Fano factor F (left panel) and the skewness C (right panel) as a function of the number 
of contacts in series N , shown for different contact transparencies. The corresponding distributions of 
transmission eigenvalues of the composite junctions are shown in fig. |5] The Fano factors Fs and the 
skewnesses Cs (normalized by e* = 2e) of the superconducting junctions are given for comparison. 
Ballistic to diffusive crossover in a disordered wire (L/l — N — 1) is shown by dash-dotted curves. 

Conclusion. We have studied the transport properties of several chaotic cavities in 
series using the circuit theory of mesoscopic transport. We obtained the analytical expression 
for the distribution of transmission eigenvalues of the composite junction as a function of 
the number of contacts and the scattering properties of a single contact. This distribution 
generalizes the previous results on noise in such a system [10] to all higher-order cumulants. 
As an example we found the first three cumulants of the charge transfer statistics both for the 
normal-state junction and in the case when one lead is superconducting. The sign of the third 
cumulant at high bias can be used to probe the asymmetry of the transmission eigenvalue 
distribution: it is negative (positive) when the weight of the distribution is more on open 
(closed) transport channels. As the number of contacts increases, all transport properties 
approach the universal limit of a diffusive wire [36] . While the crossover from a few cavities to 
the diffusive-wire limit has already been studied through the noise in the normal state [10, 11], 
experimental investigations of either higher-order correlators or cavities in contact with a 
superconductor are still to come. 
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